Fermi surface nesting and possibility of orbital ordering in FeO 
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We study FeO, a Mott insulator in GGA and GGA+U approximations. In the GGA we find a multi-band 
metallic state with remarkable inter-band nesting between two t2g bands of Fermi surface, which signals possible 
instability towards an orbital ordered insulating phase. Such broken symmetry state, although has lower energy 
than the underlying homogeneous metallic state, but the gap magnitude is less than the experimentally observed 
optical gap. Therefore we incorporate the calculated value of on-site Coulomb repulsion U on orbital ordered 
state. We find that symmetry breaking and Coulomb correlations cooperate together to stabilize the system 
and give an insulating orbital ordered state, with the gap magnitude very close to the experimental value. We 
propose this method as a possible indication of orbital ordering in LDA and GGA calculations. We check our 
method with known examples of LiV02 and LaMnOa. 
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INTRODUCTION 

Charge and spin are two fundamental characteristics of 
electrons. Pauling invented the concept of "orbital" as an- 
other important attribute of electrons. This theoretical tool 
which intuitively refers to the shape of electronic clouds, was 
so powerful that enabled a qualitative understanding of the 
electronic properties of elements in the periodic table, as well 
as binding in many molecules and solids in very early days 
of quantum mechanics, when no computers existed. Unlike 
charge and spin which directly couple to many experimental 
probes, the question of the experimental observation of orbital 
degrees of freedom remained untouched until recently 
Resonance x-ray scattering has been successfully employed 
for direct observation of the ordering of orbitals in mangan- 
ites, Ai-a^Ba^MnOs as well as some other systems such 
as transition metal oxides (See table 1 of Ref. 1). 

Observation of orbital ordering has established that orbitals 
are not just theoretical hypotheses. They are another degree of 
freedom, just like charge and spin which at low temperatures 
can form various types of ordered states. For materials such as 
KCuFa and LaMnOs, which are a priori known to be orbital 
ordered, people have employed LDA+U approximation with 
the assumption of orbital ordering to justify the ground state 
properties iH 0] • However, an important theoretical question 
is, given the Fermi surface topology and the electronic struc- 
ture of an unknown material, under what conditions the strong 
correlation effects are anticipated to stabilize an orbital or- 
dered state at low temperatures. In this paper we show that 
one can study the possibility of orbital ordering within ab ini- 
tio spin-density-functional scheme by appropriate Fermi sur- 
face analysis. We employ our method to predict the possibility 
of orbital ordered phase for FeO at ambient pressure. 

Despite the stoichiometric simplicity, the iron-monoxide, 
FeO, presents a challenge in terms of the theoretical under- 
standing of its electronic and magnetic properties iSH,!?!]. In 
series of transition metal mono-oxides the cubic crystal field 
splits the d orbitals into two-fold eg and three-fold degenerate 
t2g States. As one moves from MnO to NiO along the pe- 



riodic table of elements, the minority spins are being filled. 
In these systems, each (111) plane has majority spins aligned 
ferromagnetically within the plane (and so are the minority 
spins), while in the next (111) plane the ferromagnetic align- 
ment is in opposite direction. When we consider the filling 
of minority spins, the degeneracy left behind in t2g gives rise 
to a metallic state with partially filled t2g band for FeO and 
CoO in Density functional Theory (DFT) within generalized 
gradient approximation (GGA). However, this prediction is in 
contradiction to experiment, where it is found to be a Mott 
insulator with an optical band gap of 2.4 eV [8]. Ref. 9 sum- 
marizes the various theoretical values obtained for the optical 
gap. The Hartree-Fock approximation largely overestimates 
the gap [9] . The diffusion Monte Carlo method give the gap 
value of 2.8 ± 0.3 eV [10]. In this paper, we show that this in- 
sulating behavior can be accounted for by a orbital ordering, 
on top of which a Hubbard type of correlation further stabi- 
lizes the insulating state. The gap value we obtain is 2.2 eV. 
In this research we use QUANTUM-ESPRESSO [1 1] as DFT 
code to apply GGA and GGA+U to FeO and GGA to LiV02 
and LaMnOs. For core-electron interaction ultra-soft pseu- 
dopotentials have been used. Wavefunction and charge den- 
sity are expanded in plane- wave with 40 Ry and 400 Ry cutoff 
respectively. The Hubbard U is set to 4.3 eV for GGA+U cal- 
culation for FeO 

Nested Fermi surfaces can potentially lower the symmetry 
of the ground state by either developing some sort of long- 
range order [12], or appropriate distortion of the lattice struc- 
ture 1I13I1 . Efremov and Khomskii 1I14I1 developed a theory 
of orbital ordering in terms of inter-band nesting between the 
two Cg bands in manganites. According to their work, the un- 
derlying inter-band nesting can lead to instability towards a 
symmetry broken phase with ordering in Cg orbitals. In this 
work we analyze the GGA Fermi surface of FeO. We observe 
a remarkable nesting pattern between two (9/ the three it2g) 
bands crossing the Fermi level. This qualitative observation 
can be made quantitative by calculating appropriate inter-band 
susceptibility which peaks around q = 1/3G, where G is 
the reciprocal lattice vector in (111) plane. The evidence for 
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FIG. 1: (Color online) Fermi surface for (a) FeO, (b) LiV02 and (c) LaMnOs- The Fermi surfaces are plotted in the reciprocal lattice (not in 
the first Brillouin zone). For FeO and LiV02 the lattice is trigonal and for LaMnOs is tetragonal. 



instability towards orbital ordered phase from GGA metallic 
ground state of FeO is substantiated with another independent 
GGA calculation for the broken symmetry phase. We also 
perform a GGA+U calculation to compare states with bro- 
ken symmetry and those without broken symmetry. Thus, or- 
bital ordering (00), and Hubbard correlations work together 
to give a Mott insulator with orbital ordering pattern, in agree- 
ment with a similar scenario proposed for LiV02 1 15]. 

To verify our approach, we apply this analysis to the 
well known orbital ordered materials, such as, LiV02 and 
LaMnOs. The former example is closer to FeO in two re- 
spects: (i) The ordering pattern occurs in a triangular lattice 
of (111) plane, (ii) Both LiV02 and FeO have degenerate 
t2g bands; however, the difference is that in LiV02 we have 
a dp class fTs], while the minority spin electrons of FeO be- 
long to class. Note that the majority spin states are well 
separated from those of minority spins |7]. In case of FeO 
we are dealing with Fe^+ cations with 3(i^ configuration. By 
Hund's rule, such configuration is equivalent to d^ class in 
the minority spin sector. We first calculate the GGA bands 
and Fermi surface for the trigonal structure. This gives us 
a metallic state with remarkable inter-band nesting shown in 
Fig- Ea). In this figure, the nesting is only between two of 
the t2g Fermi surfaces. The third t2g band does not contribute 
substantially to the Fermi surface. We have checked that this 
band is further pushed away from the Fermi level by small 
amounts of on-site Coulomb interaction U. Hence despite the 
threefold degeneracy of t2g bands, here we are dealing with a 
two-band situation to which the picture proposed by Efremov 
and Khomskii in the context of orbital ordering instability of 
Cn bands can be applied. The same type of inter-band 
nesting can be observed in the GGA Fermi surfaces of LiV02 
(Fig-IHt))), which is known to exhibit orbital ordering at low 
temperatures lE5l[l3,[i3]- The example of LiV02 is very sim- 
ilar to FeO, in terms of underlying triangular lattice structure 
in (111) plane. This is also reflected in their Fermi surface 
topologies in Fig.Q In Fig.IHc) we plot the GGA Fermi sur- 
face of LaMnOs. As shown by Efremov and Khomskii [14], 
such a nesting leads to an instability towards orbital ordered 



state at low enough temperatures. In the following we are go- 
ing to argue that the same argument predicts orbital ordering 
for LiV02 and FeO. 

To quantify the above discussion, we calculate the contribu- 
tion of inter-band processes to generalized susceptibility de- 
fined by. 



Xinter(q) = ^ 



/(^£)[i-/(4+J] 



(1) 



where a and h denotes the two nested bands, /(e) is the occu- 
pation number and e^'^ denotes the band energy. To identify 
the exact value of nesting vector we calculated generalized 
susceptibility in direction of Gi ± G2, where Gi and G2 are 
two basis vectors in the reciprocal space. 




FIG. 2: (Color online) Generalized susceptibility for FeO, LiV02 
and LaMnOs. For FeO and LiV02 the direction of q is parallel 
to Gi — G2 where Gi and G2 are vectors in reciprocal space of 
the trigonal lattice. For LaMnOs we chose q in the direction of 
Gi + G2 where Gi and G2 are reciprocal lattice vectors of the 
tetragonal lattice (| G-^ |= | G2 |). The vertical line indicates the 
maximum of generalized susceptibility at q = | (Gi — G2) for FeO 
and LiV02. 

In Fig. [2] we plot the results of numerical evaluation of inte- 
gral ([B for values of q in the direction indicated by arrows on 
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FIG. 3: (Color online) Charge density of minority spin of Fe atoms 
in (111) surface with GGA calculation. The iron atoms are slightly 
displaced and then they are allowed to be relaxed. As can be seen, the 
new minimum in this GGA calculation corresponds to an orbital or- 
dered state. The new cell shown in this figure is y/Sa x VSa inspired 
by the nesting vector obtained from the generalized susceptibility. 




FIG. 4: (Color online) Orbital ordering in the GGA+U solution. This 
pattern is analogue of the mean field state proposed in Ref.FTsI 
for LiV02. As can be seen the incorporation of correlation effects 
via [/, does not destroy the orbital ordering. This figure corresponds 
to the state with (a, (3, 7) = (1, 0, 0). 



panels (a)-(c) of Fig.[TJ respectively [18]. As can be seen in 
Fig.O for the case of LaMnOa, the inter-band susceptibility 
peaks around the ordering vector, q = ^(Gi + G2) in (001) 
plane |14]. This implies an ordering in the xy plane with 
2a X 2a unit cell. Now let us focus on the case of t2g systems, 
FeO and LiV02. As can be seen in Fig.[2l for the case of FeO 
and LiV02, the susceptibility peaks around q = |(Gi — G2). 
Based on the theory of Ref. lll4l1 . we speculate this implies a 
possible instability towards ordering of t2g bands on a trian- 
gular lattice of (111) plane. 

Feeding these information back into the DFT machinery, 
we examine the total energy, gap magnitude and charge dis- 
tribution profile assuming an ordering pattern suggested from 
our Fermi surface analysis. Hence we are lead to examine a 
broken symmetry phase with a >/3a x \/3a supercell in (111) 
plane The structure contains three Fe atom per (111) 

plane. With the new supercell, we allow for a small displace- 
ment of Fe atoms (according to Ref. 7) to check for the pos- 
sibility of relaxation towards a new minimum in the GGA ap- 
proximation. This new state is orbital ordered insulating state, 
as can be seen in the charge density profile shown in Fig. [3l 
The insulating behavior of this state is due to broken symme- 
try accompanying the orbital ordering. The gap magnitude at 
F point for this state is E^^^ :^ 0.54 eV. The energy of this 
new minimum per supercell is about 0.47 eV lower than the 
corresponding non-ordered metallic state. 

Although in the GGA approximation when we allow for 
the possibility of a broken symmetry state, we find an insu- 
lating state with lower energy than the original metallic one, 
but the energy gap is still far less than the experimentally ob- 
served value of = 2.4 eV LSJ. The GGA energy level of 
FeO suggests the following picture: Five majority spins are 
well below the Fermi level, leaving one minority spin in the 
relevant t2g band IlZJ . To this extent, we consider FeO to be- 
long io class. Therefore we are dealing with one minority 
spin per site Mott insulating situation. Hence, we have to take 
into account the role of on-site Coulomb interaction, U. We 



incorporate the effect ofU ~ 4.3 eV within the GGA+U ap- 
proximation. With this value of [/, we obtain the energy gap 
^GGA+u _ 2.2 eV at the F point. Note that we did not ad- 
just the value of U. The above value for U was calculated 
in Ref. 7. Cornpared to sophisticated quantum Monte Carlo 
calculations lldl , our result is in remarkable agreement with 
the experimental value fs^. Therefore both orbital ordering, as 
well as Coulomb correlations are important to understand the 
electronic structure of FeO Mott insulator. 

With GGA+U there is a freedom on choosing the appro- 
priate state in the t2g subspace to be occupied by minor- 
ity spin. The most general state in this subspace is of the 
from, a\xy) + I3\yz) + 7|zx). The underlying triangular 
lattice structure with C3 symmetry requires invariance un- 
der 120° rotations. Therefore the two other states are re- 
quired to be, a\yz) + I3\zx) + ^\xy), a\zx) + I3\xy) + ^\yz). 
Each set of values for (a,/3, 7) corresponds to a new or- 
bital ordered configuration. The only choice which leads to 
state without orbital ordering is the homogeneous one with 
W\ = \P\ = ItI = 1/^3 (the so called aig state). In GGA+U 
it is not feasible to check the energy of all possible states. In 
table II we compare the energy of the non-ordered aig state 
with two states corresponding to (a,/3, 7) = (1,0,0), and 
(a,/3,7) = (1,1,0)/ V2. Second column of table [B indicates 
that, in presence of electron-electron interaction U, orbital or- 
dering still reduces the total energy by tens of meV per unit 



TABLE I: Energy gained by orbital ordering with respect to non- 
ordered aig state. For each case we include the energy of system 
without atomic position relaxation and with atomic position relax- 
ation. The energies are expressed as energy per unit cell (each cell 
contains six Iron atoms). 



Configuration 


E-Ea,g(meV) 


E-Eaig(meV) 


(a,/3,7) 


without relaxation 


with relaxation 


(1,0,0) 


-60 


-174 




-36 


-66 
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cell. The state (1, 0, 0) has lower energy for which we have 
plotted the charge density in Fig.lH This state is the negation 
of the state proposed in Ref. 1 15] for the LiV02. Such a nega- 
tion image is reasonable given the fact that and d'^ classes 
in t2g subspace are connected with a particle-hole symmetry 
transformation. 

So far we have checked that electronic interactions can take 
advantage of the instability suggested by the nested Fermi sur- 
face, to stabilize an orbital ordered pattern. Next we ask the 
question, can electron-lattice interactions take advantage of 
this instability to give rise to appropriate form of Jahn-Teller 
distortion? To verify this for the case of FeO, we start with 
orbital ordered states (a, 7), then we allow the atomic po- 
sitions to relax. In the third column of table ID we report the 
calculated total energy per unit cell when the relaxation is al- 
lowed. As it can be seen, the relaxed structure has lower en- 
ergy than the corresponding non-relaxed state. Note that to 
verify whether the orbital ordered state is more stable than 
the non-ordered state, we first distorted the atomic positions 
which resulted in orbital ordering pattern of Fig. O But here 
we check the reverse sequence, i.e. we start off with an or- 
bital ordered state, then we check whether the atomic dis- 
placements can further stabilize orbital ordered state or not. 
Therefore for both orbital ordered states considered here, the 
ordering is always accompanied with cooperative Jahn-Teller 
distortion. 

To summarize, we started with a GGA metallic band pic- 
ture. The inter-band nesting pointed us to examine the en- 
ergy of orbitally ordered state within the same approxima- 
tion. However, since in the case of transition metal oxides 
the Coulomb correlation U is also important, we also com- 
pared the energy of orbital ordered states in presence of on- 
site Coulomb interaction U, in the GGA+U approximation. 
In FeO, Orbital ordering significantly stabilizes both GGA 
as well as GGA+U states with respect to corresponding non- 
ordered phases. We also allowed for relaxation of the atomic 
positions on top of orbitally ordered GGA+U state. We found 
that the Jahn-Teller distortion accompanying the orbital order- 
ing is automatically realized in our approach. As a check of 
this method, we also studied the cases of LiV02 and LaMnOs 
which are known to exhibit orbital ordering at low tempera- 



tures. Therefore we propose this Fermi surface analysis as a 
rout to explore the possibility of orbital ordering and/or re- 
duction in the lattice symmetry, within the DFT (LDA/GGA) 
electronic structure calculations. 
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